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Abstract. We obtain various new well-posedness results for continuity and transport equa- 
tions, among them an existence and uniqueness theorem (in the class of strongly continuous 
solutions) in the case of nearly incompressible vector fields, possibly having a blow-up of the 
BV norm at the initial time. We apply these results (valid in any space dimension) to the 
k X k chromatography system of conservation laws and to the k x k Keyfitz and Kranzer 
system, both in one space dimension. 



1. Introduction 

1.1. Continuity and transport equations. In the last few years some progress has been 
made on the well-posedness of the continuity and transport equations dtp + div (bp) = 0, 
dfW + b ■ Vw = 0, for weakly differentiable (in the space variables) velocity fields b. The 
first seminal paper by DiPerna and Lions [22] considered Sobolev vector fields with spatial 
divergence div 6 in L^([0, T]; L°°), and more recently in [1] Ambrosio extended the result to 
BV vector fields, assuming absolute continuity of the divergence and -^^^([0, T]; L°°) regularity 
of its negative part (see also the lecture notes by Ambrosio and Crippa [1]). This bound 
is the simplest way, as it can be easily seen with the method of characteristics, to prevent 
blow-up in finite time of the solutions; on the other hand, it was already clear from the 
application made in 0, [3] to the Keyfitz and Kranzer system [2H] , that this assumption is 
not very natural. The Keyfitz and Kranzer system has the peculiar form 

dtu + div {uf{\u\))=0 
and (at least formally) it decouples in a scalar conservation law for p := namely 

a,p + div(p/(p)) = 0, 
and a transport equation for the "angular" part 6, related to m by -u = p6: 

dte + f{p)-ve = o. 

If we take p as the entropy admissible solution of the scalar conservation law, it turns out that 
the velocity field b = f{p) in the transport equation belongs to BV if the initial condition 
p is in BV, but its distributional divergence needs not be absolutely continuous; the only 
(weaker) information available is that there is a bounded function, namely p, transported 
by b. This leads to the concept of nearly incompressible vector fields (namely those b for 
which a bounded transported density p exists), and well-posedness results in this class of 



2 



LUIGI AMBROSIO, GIANLUCA CRIPPA, ALESSIO FIGALLL AND LAURA V. SPINOLO 



vector fields have been investigated in [6], [181 120], also in connection with a compactness 
conjecture made by Bressan [13]. See also [T7j for a counterexample to the applicability of 
these techniques to general multidimensional systems of conservation laws. 

In this paper we obtain a well-posedness result for this class of vector fields, see The- 
orem 12.51 for a precise statement. However, the main result is an improvement of the 
L^{[0,T]; BVioc) condition considered in all previous papers: going back to the case of the 
Keyfitz and Kranzer system, it is clear that we cannot expect this regularity for p (and then 
for b = f{p)), unless p G BV. More precisely, in the one-dimensional case and imposing 
suitable conditions on the flux function /, if p G L°° we gain indeed a regularizing effect 
(see Oleinik [22]). However, this provides only L^^^(]0,T]; BVioc) regularity, and also (using 
the equation satisfied by p) BVioc{]0,T] X R'^) re gularity. This last assumption on the vector 
field is critical in view of a non-uniqueness example provided by Depauw [2T], where the 
vector field has precisely this regularity. Our main result, given in Theorem 12.71 is that for 
this class of vector fields (adding bounds on the divergence in the same spirit of the nearly 
incompressibility condition, see Definition 12. 6p existence and uniqueness can be restored, 
provided one works with strongly continuous in time solutions. Here, of course, the main 
difficulty is in the existence part, since standard approximation schemes in general do not 
provide this strong continuity property. 

The main application of our results concerns the so-called chromatography system. As 
discussed in Section II. 2^ we obtain new well-posedness theorems for this equation. As a 
byproduct of our analysis, we also have applications to the one-dimensional Keyfitz and 
Kranzer system, obtaining in particular a well-posedness theorem for bounded initial data 
(results in this flavor were already known: see for example Freistiihler [23] and Panov |30j). 

1.2. Chromatography. Although all our discussion and results can be easily extended to 
the k X k chromatography system (see Remark 14.101) . for simplicity of exposition we will 
always consider only the case k = 2. The chromatography equation 



is a system of conservation laws belonging to the so-called Temple class. Consider first a 
general system of conservation laws 



The associated Cauchy problem is well-posed under the assumption that the total variation 
of the initial datum U is sufficiently small (see e.g. [12] and the references therein). If the 
total variation of U is bounded but large, then the solution may experience blow up in finite 
time, as shown for instance by Jenssen [25] . 






dtU + d,F{U) = 



where {t,x) G [0,+oo[xM and U G M''. 
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Temple systems were introduced in |33] and they are defined by special properties imposed 
on the structure of the eigenvector fields of the Jacobian matrix DF{U). Thanks to these 
features, well-posedness results for Temple systems are available for a much larger class 
of initial data compared to general systems of conservation laws. In particular, Serre [3l] 
obtained global existence of weak solutions for 2x2 Temple systems with no smallness 
assumptions on the total variation of the initial datum. By relying on wave front-tracking 
techniques, in [8] Baiti and Bressan constructed a Lipschitz continuous semigroup defined on 
initial data with large but bounded total variation. The extension of the semigroup to L°° 
initial data was first achieved by Bressan and Goatin [H] under a nonlinearity assumption on 
the structure of the eigenvector fields of the Jacobian matrix DF{U), and then extended by 
Bianchini in [9] to general Temple class systems, by removing any nonlinearity assumption. 
One of the main difficulties faced in [5] is that if the nonlinearity assumption fails, one 
cannot hope for a Lipschitz dependence of the semigroup on the initial data. This is shown 
by the example discussed in [T3], which involves precisely a family of Cauchy problems for 
the chromatography system. 

The main result in [9] is the existence of a semigroup St of weak solutions defined on L°° 
initial data and satisfying the following conditions: 

(1) (t, U) I—*- StU is continuous with respect to M x topology, with values in L^^^; 

(2) If U is piecewise constant, then for t sufficiently small SfU coincides with the function 
obtained by gluing together the "standard" solutions constructed as in Lax 

The proof exploits the so called wave front-tracking algorithm, which relies on the approx- 
imation of the initial datum with piecewise constant functions. In [9] uniqueness is also 
obtained: namely, it is shown that there exists a unique semigroup of weak solutions satis- 
fying conditions (1) and (2) above. In addition, the weak solutions provided by the maps 
t I— StU are automatically entropy admissible. 

The results in [9] apply to general systems in the Temple class, while in the present work 
we restrict to the chromatography system. The main novelties here are the following. 

First, our approach is completely different from the one in [H]. Namely, by introducing a 
change of variables in analogy to Ambrosio and De Lellis [S] and to Bressan and Shen |15j . 
we split the chromatography system in the coupling between a scalar conservation law and 
a transport equation, and then we heavily exploit transport equation techniques. 

The second point concerns uniqueness: we manage to get it in the classes of functions 
defined in SectionlHwithout requiring any stability with respect to perturbations in the initial 
data. Conversely, the approach in [9] is in the spirit of the Standard Riemann Semigroup (see 
Bressan [Uj): one shows that there exists a unique semigroup of weak solutions satisfying 
(1) and (2) above, thus one needs to prove stability and the "right" behavior on piecewise 
constant data to get uniqueness. This does not correspond, a priori, to requiring uniqueness 
of weak entropy admissible solutions. Finally, we point out in passing that we manage to 
relax the hypothesis of strict hyperbolicity made in [9], which fails when the initial data Ui 
and U2 vanish in some region (see Definition 14. 3p . However, the price we have to pay in this 
case is the requirement -ui + ^2 G BVioc(JR)- 
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The idea of attacking the chromatography system via a change of variables which sphts 
the system into a couphng between a scalar conservation law and a transport equation 
came from Bressan and Shen [15] (see also Panov [301 Remark 4]). In [15] the authors are 
mainly concerned with the study of ODEs with discontinuous vector fields: they formulate 
conditions under which uniqueness holds, and they show that the vector field obtained from 
the chromatography system after their change of variables satisfies these conditions, provided 
the initial data take value in a particular range. In our proof we perform a different change 
of variables, which has the advantage of being linear, and hence behaves well under weak 
convergences. In particular, we can come back to the original system, obtaining distributional 
solutions which are admissible in the sense described in the Appendix. Also, we manage to 
handle a larger class of initial data. 

Acknowledgement. We thank Fabio Ancona and an anonymous referee for pointing out 
to us the extension of our results to k x k systems, as illustrated in Remark I4.10[ 

1.3. Content of the paper. In Section [2] we prove the basic well-posedness results for the 
continuity equation, first in a class of bounded nearly incompressible vector fields having 
-BVioc([0,T] X Mf^) regularity, and then under additional assumptions on the transported 
density p (see Definition 12.61) for bounded vector fields having BVioc(]0,T] x M'^) regularity. 
In Section [3] we discuss Depauw's example [21] and its relation with our well-posedness 
results. In Section H] we present an application of these results to the chromathography 
system. First we classify in Lemma 14.11 the entropy-entropy flux pairs for this system, and 
then we prove two basic existence and uniqueness results: the first one, Theorem l4.4t provides 
existence and uniqueness in L°° under the assumption that the sum ui + U2 of the initial 
conditions is in BVioc] the second one. Theorem 14. 8^ replaces this regularity condition with 
infj<"Ui + ^2 > for any K G M. compact. The applications to the Keyfitz and Kranzer 
system are discussed in Section [5j in Theorem 15.21 we show existence and uniqueness of 
strongly continuous renormalized entropy solutions, when the initial datum u is bounded 
and satisfies infi^ \u\ > for any K cM. compact. In the Appendix we list a few basic facts 
on scalar conservation laws needed in the paper. 

1.4. Main notation and conventions. Finally, we specify some (standard) convention 
about spaces of time-dependent functions: if J C M is an interval and X is a separable 
Banach space, by U'{J]X) (resp. L\^^{J]X)) we mean all measurable functions u : J ^ X 
such that \\u\\x G L^iJ) (resp. L^il) for all / C J compact). In the particular cases of 
non-separable spaces X, in this paper X = L°° or X = BV, the definition is the same, but 
measurability is understood in a weak sense, using the embedding of these spaces in Ll^^. 

We shall often consider locally bounded (in space and time) distributional solutions to 
the continuity equation dtp + div (bp) = 0, with b locally integrable. In this case it is well 
known (see for example Lemma 1.3.3 in [TH]) that the map t i— > p{t, ■) has a unique locally 
weakly-* continuous representative (i.e. t ^ j p{t,x)(j){x) dx is continuous for any G L°° 
with compact support) and we shall always work with this representative, improving in some 
cases the continuity from weak to strong. We shall use the notation L°° — w*, L} — s for the 
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weak-* and strong topologies, and L"^^ — w*, L/^^ — s for their local counterparts. Finally, 
we will often consider the set of functions which are continuous in time with values either in 
L°°(M'^) — w* or in L^(M'^) — s (or in their local counterparts). These spaces are denoted by 
C°([0, +oo[; L°°(M'^) - w*) and C°([0, +oo[; L^M'^) - s), respectively. 

2. WELL-POSEDNESS results for CONTINUITY EQUATIONS 

In this section we study the well-posedness of the continuity equation 

dtu{t, x) + div [b{t, x)u{t, x)) =0 
m(0, x) = u{x) . 

Our main results are presented in Theorems 12.51 and 12.71 

2.1. A preliminary renormalization lemma. We first prove a technical lemma of stan- 
dard flavour, regarding renormalization and strong continuity for solutions u of equation 
(12. 2p with a BV nearly incompressible vector field, see De Lellis [20] for a systematic treat- 
ment of this topic. The main difference is that, for our subsequent discussion, we need to 
consider functions q which are possibly not bounded away from zero, or even equal to zero 
on a non-negligible set. 

Lemma 2.1. Let b G i?Vioc([0, +oo[xM'^; M'^) be a bounded vector field. Assume that there 
exists a nonnegative locally bounded function 

q G 5Moc([0,+oo[xM'^) nC°([0,+oo[; L^i^'') - w*) 

satisfying \Dq{t, ■)\{Bj^) G L^dl'^, +'^[) for all R> and 

dtq + div (bq) = 

in the sense of distributions in ]0, +oo[x]R'^. Let u G LJ^^([0, +oo[x]R'^) be a distributional 
solution of 

dt{qu) + div [bqu] = 0. (2.2) 
such that the map t qu{t, ■) is weakly-* continuous in L^^{M.'^). Then: 

(i) u is a renormalized solution, in the sense that the map 

t^q{t,-)P{u{t,-)) 
provides a distributional solution of 

dt{qf3{u))+ div {bq(3{u)) =0 (2.3) 

on ]0, +oo[x]R'^, for any function [3 G LipiQj,(M). 

(ii) // in addition the map 

t^q{t,-) 

is strongly continuous from [0, +oo[ with values in Ll^^lM."^), then also the map 

t (-^ qu(t, ■) 

is strongly continuous from [0, +oo[ with values in Ll^^{M.'^). 



6 LUIGI AMBROSIO, GIANLUCA CRIPPA, ALESSIO FIGALLL AND LAURA V. SPINOLO 

(iii) Assume that ui and U2 G Lj^j,([0, +oo[xM'^) are solutions of (\2.2\i such that there 
exists r > satisfying 

g(r, x)Mi(r, x) = q{T, x)u2{t, x) for a.e. x E . 

Then for all t G [0, +oo[ we have 

q{t,x)ui{t,x) = q{t, x)u2{t, x) for a.e. x G W^. 

Proof Let us define B{t,x) := {q{t,x),b{t,x)q{t,x)). Then B G Wioc([0, +oo[xM^; M x M^) 
and is divergence free in space-time. We introduce an artificial "time variable" s G [0, +oo[ 
and observe that the function 

v{s, t, x) := u{t, x) 

is a distributional solution of 

dsv + divt,^. (Bv) = on ]0, +oo[,x]0, +oo[txM^. 

Since B{t,x) is an autonomous (i.e. independent of s) divergence free vector field belonging 
to L°°([0, -|-oo[s; i?Vioc([0, +oo[iXM^; MxM'^)) , thanks to Ambrosio's renormalization theorem 
(see [1], or [21 Section 5]) the function P{v{s,t,x)) solves 

dsp{v) + divi,^ (Bpiv)) =0 on ]0, +oo[x]0, +oo[xM'^, 

or equivalently 

dt{qP{u)) + div {bqp{u)) = on ]0, +cx)[xM'^. 

This proves (i). 

We now show (ii). Applying the result in (i) with (3{u) = m^, we have that q{t, ■)[u(t, ■))^ 
provides a distributional solution of (12. 3p . Hence, there exists z G C°([0, +oo[; L'^^(W^) —w*) 
such that 

z{t,-) = q{t,-){u{t,-)Y a.e. in M*^, 

but only for almost every t > 0. Since t (l{t,-) is strongly continuous with values in 
Li^^(M"'), it follows that 

g(t, ■)z{t, ■) G C° ([0, +oo[; L^M') " • 

Notice that the difficulty in the proof of the strong continuity oft ^ q(t, ■)u(t, ■) comes from 
the fact that it is not a priori obvious that the equality 

q{t,-)z{t,-) = {q{t,-)u{t,-)Y , 

which is true for almost every t > 0, in fact holds for every t > 0. This will be shown by 
truncating the equation at an arbitrary time. 

Let us fix to G [0,+oo[: our goal is showing that t t-^ qu{t,-) is strongly continuous in 
Lj'Q^(M"') at to- Let us define the functions 

' ^ \ b{t,x) if t G]to,+oo[, 
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q{to,x) ifte]-oo,to] 
q{t,x) if t e]to,+oo[, 

u{tQ, x) if t g] — oo, to] 
u{t, x) if t E]to, +oo[ 

and notice that q(to, x) G BViodJ^'^) by the weak continuity of t i— > q(t, ■) and the uniform in 
time and local in space bound on \Dq{t, ■)\. Thanks to the weak continuity of q and qu we 

get 

dtqto + div (btoqto) = and dt{qtaUto) + div [h^qtoUt^,) = 

in the sense of distributions on M x M*^. By applying the same strategy as in the proof of (i) 
with the vector field Btg := {qto , h^qto) , we deduce that 

dt{qt,ul) + div {bt,qt,ul) =0 on M x M"' 

in the sense of distributions. Combining this with the fact that the map t qto{t, ■)uf^{t, ■) 
is constant and equal to (gn^)(to, ■) for t < to, we easily deduce that 

qu'^{t, ■) qu^ito, ■) 

weakly-* in Lj^^(]R'^) ast^t^, and thanks to the strong continuity of t i-^ ■) this implies 

gV(t,-)^gV(to,-) (2.4) 

weakly-* in L;^^(R'^) ast^f^. 

From (12.41) combined with the weak-* continuity in Lj^j,(M'^) of the map t qu{t, ■), we 
deduce that 

quit,-) qu{to,-) 

strongly in Ll^^iM.'^) as t ^ t^. This proves the right continuity at to. The proof of the left 
continuity is analogous. 

Finally, to show (iii), we observe that vi{s,t,x) := ui(t,x) and V2{s,t,x) := U2{t,x) are 
both solutions of 

dsv + divt,^. {Bv) = on ]0, +oo[,x]0, +oo[txMf.. 
Hence, by Ambrosio's renormalization theorem, \vi — solves 



ds{\vi -V2\) +divt,^. [B\vi -V2\) =0 on [0, +oo[sX [0, oo[ 



or equivalently 

dt{q\ui - U2\) + div [bq\ui - U2I) = 0. (2.5) 

Moreover, by considering the functions bt^^, qt^ and {u2 — ui)to as in the proof of (ii), we can 
apply as above the renormalization with f3{u) = \u\ to obtain 

dt{qto\{ui - U2)to\) +div {btoqto\{ui - U2)to\) =0, 
which as above implies that the map 

t ^ q{t, ■)\u^{t, ■) - U2{t, ■)\ G L~ (M'^) (2.6) 
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is weakly-* continuous. Hence, arguing as in [201 Lemma 3.17], (12. 5p and (12.61) imply 

/ q{t, x)\ui{t, x) — U2{t, x)\ dx < / q{s, x)\ui{s, x) — U2{s, x)\ dx 

for all t, s > 0. In particular, by setting s = r, we obtain the desired result. □ 

Remark 2.2. Assume that in Lemma [2.11 we replace the hypothesis of strong continuity of 
1 1— *• g(t, ■) with the condition that for every t 

^ < q{t,x) < C a.e. x G M"^ 

for some constant C > 0. Then the map t i— * u{t, ■) is strongly continuous in Ll^^. This is 
shown in De Lellis [201 Corollary 3.14]. 

2.2. Well-posedness of weakly continuous solutions. In this section we prove a well- 
posedness result for (12.11) in a function space adapted to the vector field b: this space depends 
on an auxiliary function p linked to b by (12. 7p . 

Definition 2.3. Let f2 be an open set in an Euclidean space and let p G L^^(jl) be non- 
negative. We denote by L°°{p) the set consisting of the measurable functions w defined on 
Q such that there exists a constant C satisfying 

|wi(y)| < Cp{y) for a.e. y E fl. 
The smallest constant C will be denoted by ||u;||j;^oo(p). 

Remark 2.4. Notice that, for any w G L°°{p), we necessarily have w = a.e. on {p = 0}. 
Hence, by adopting the convention that = if p{y) = 0, the quotient 

p{y) 

is well-defined on Q and is a bounded measurable function satisfying (w/p) ■ p = w a.e. 

Theorem 2.5. Let b e BVioc{[0, +oo[xW^;W^) be a bounded vector field. Assume that there 
exists a nonnegative locally bounded function p G i?Vioc([0, -|-cx)[xM'^) satisfying |-Dp(t, ■)\{Br) G 
Li^^([0, +oo[) for allR>0 and 

dtp + div {bp) = (2.7) 

in the sense of distributions. Then: 

(i) For any initial datum u G L°°{p{0, ■)) there exists a solution u G ^^([0, +oo[; L°°(M"') — 
w*) nL°°{p) to the Cauchy problem (12.11) . 

(ii) The solution u to (12. ip is unique in the class C°([0, +oo[; L°°(M^) - w*) n L~(j9). 

(iii) \\u{t, ■)\\L--{p{t,-)) < ||^^||l°°(p(o,-)) M all t > 0. 
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Proof. Properties (i), (ii) and (iii) are trivially satisfied when p is identically zero. Hence 
in the following, without any loss of generality, we assume that p is not identically zero. 
Existence. Let us define /q := u/p{0, ■) (see Remark [2. 4p . Roughly speaking, the strategy 
to construct a solution is to solve 

dtl + b-Vl = 
/(0,-) = /o, 

and to define u = Ip. Since b is not smooth, we need a regularization argument. 

We first remark that the assumption that p is nonnegative and not identically zero implies 
that p{t, ■) is not identically zero for any t > 0. Indeed, assume by contradiction that 

p(r, x) = for a.e. x G M*^ 

for some r > 0. Then by applying Lemma [2.1 ( iii) with q = p, ui = 1 and M2 = we obtain 
that for every t > 

p{t,x) = for a.e. x G M'^, 

against our assumption. 

We now consider a sequence 77^ of convolution kernels supported on the whole M*^ (Gaussian 
kernels, for instance), and we define 

, (bp) * Ve 

Pe.= p*r]e, be'.^ 



P*Ve 



Then p^ is smooth and strictly positive everywhere, so that b^ is smooth and bounded. We 
now solve for every e > the Cauchy problem 

dtls + be-Vle = 



Then, using the identity dtPe + div (bePe) = 0, one can easily check that 

r dt{hPe)+diY{bJ,p,)=0 
\ lePe{0,-) = U*ri^. 

Since 

IKe(^, ■)||l°°(M'') < \\h{^,-)\\L^(^d) < ||^i||L-(p(o,-)) for every t G [0,+oo[, (2.8) 



the functions le are uniformly bounded in ([0, +oo[x]R'^) . Hence, up to subsequences, 
converges weakly-* in L°°([0, +(X)[xR"') to a function / G L°°([0, +oo[xR'^). Observing that 
Pe ^ p in l^ioc([0, +00 [xM*^), we get that 

lePe ^lp:=u weakly-* in L'^^{[Q, +oo[xM'^) , 
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and u solves (12.11) . Moreover, thanks to (12. 8p . (iii) holds. Finally, the fact that /^Pe are 
uniformly continuous in time with respect to the weak-* topology of L5^^(M'^) (see for in- 
stance [71 Lemma 8.1.2]) implies that u belongs to C°([0, +oo[; L°^(]R^) —w*^ and u(0, ■) = u. 
This proves (i). 

Uniqueness. Let ui, U2 G L°°{p) be two solutions of (12.11) . We apply Lemma [211 111) 
with q = p and r = to the functions Ui/p and U2/P (remember Remark 12. 4p . obtaining 
that, for all t G [0, +00 [, 

p(t, x) ^^4^ = Pit, ^^4^ for a.e. X e M". 
p{t,x) p{t,x) 

Recalling again Remark [231 this gives that for all t G [0, +00 [ the equality Ui(t, x) = U2(t, x) 
holds for a.e. x G M"^, proving (ii). □ 

2.3. Well-posedness of strongly continuous solutions. In this section we weaken the 
assumption b G 5Vioc([0, +oo[xM'^; M*^) replacing it by 6 G 5Vioc(]0, +oo[xM'^; M'^) and, 
by adding some conditions on the density p transported by b as in (12.71) . we can still get 
well-posedness of (12. ip in the class of strongly continuous solutions. 

Definition 2.6. Let b : [0, +oo[xM'^ — M'^ be a bounded vector field. We say that b is locally 
strongly nearly incompressible if there exists a function p : [0,+oo[x]R'^ R such that the 
following properties are satisfied: 

(1) t (-+ p{t, ■) is strongly continuous with values in Ll^^(M.'^); 

(2) For every R> and T > there exists a constant Cr^t > such that 

^ < pit,x) < Cr^t for a.e. (t,x) G [0,T] x 5^(0); (2.9) 



(3) The equation 

dtp{t, x) + div {b{t, x)p{t, x)) = (2.10) 
holds in the sense of distributions on ]0, +00 [xM'^. 

Theorem 2.7. Let b G i?Vioc(]0, +oo[x]R'^; M'^) satisfy the assumptions of Definition \2.6[ for 
some function p which in addition belongs to BVioc{]0, +00 [xM'^) and satisfies \ Dp{t, ■)\{Br) G 
-^ioc(]0) +^[) foi^ all R > 0. Then there exists a locally bounded solutionu G C°([0, +oo[; i^ioc(K 
s) to the Cauchy problem (12. ip . Furthermore the solution is unique in this class. 

Proof. Uniqueness. Let mi, U2 G C°([0, +oo[] L\^^{W'-) - s) be locally bounded solutions of 
(12. ip . We first observe that, from the strong continuity, it follows that 



/ \ui{e,x) — U2{e,x)\dx ^ Q as e — O"*" (2.11) 



for every > 0. 
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We now apply Lemma [2 .1( 1) with i3{z) = \z\ and q = p to the function u = {ui — U2)/ p, 
on the time interval [r, +00 [, for an arbitrary r > 0, to deduce that 

dt[\ui - M2I) + div [b\ui - M2I) = (2.12) 

in the sense of distributions in ]0, +oo[xM'^. Arguing as in [201 Lemma 3.17], by fl2.12p we 
easily obtain that for any 0<£:<t<T<oo and for any i? > 0, there holds 



Ui{t,x) — U2{t^x)\dx < I \ui{e,x) — U2{s,x)\dx . (2-13) 

By letting e ^ 0^ and recalling (12. lip , thanks to ( 12.13^ we get that, for any t G [0, +oo[, 

x) = M2(t, x) for a.e. a; G -B_r(0). 

Since R is arbitrary, we conclude that ui = U2- 

Existence. The proof is organized in three steps. 

Step 1. We describe the issues we have to address and we sketch how we will proceed in 
the remaining two steps. 

Under the assumptions of the theorem, we can apply Lemma [2?lT ii) to the solution u/p on 
the domain [r, +00 [xM^ for any r > and deduce that, for every locally bounded solution 
u of fl2.ip . the map 

1 1— > u(t, ■) 

is strongly continuous from ]0, +00 [ with values in LJq(,(M'^). Thus the only issue is construct- 
ing a solution which is also strongly continuous at t = 0. 

The strategy is first approximating in Ll^^ our vector field b with smooth strongly nearly 
incompressible vector fields be with corresponding smooth densities p^. Then we observe 
that, if Zs is a solution to the continuity equation 

dtZe + div (beZs) = 
Ze{0, ■) = n * , 

then Ag := ^^g/pe is a solution to the transport equation 

dtle + be ■ VA, = 

A.(0,-)-^. 

Noticing that Pe(t, ■) is strongly continuous by assumption, we only need to prove the conti- 
nuity at of t t-^- ■). This can be directly shown by using the representation formula for 
the solution and exploiting the local strongly nearly incompressibility (see Definition 12. 6p . 
Finally, we exploit the weak-* compactness of z^, we let e — > 0, and we prove that any limit 
point of provides a strongly continuous solution. 

Another issue we need to address is that the vector field b is only locally nearly incom- 
pressible, meaning that the bounds (12. 9p are not necessarily global. We thus introduce a 
localization argument, by first restricting to initial data with compact support: because of 
the finite propagation speed, the solution we construct has then compact support for every 
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t > 0. We finally obtain the solution for general initial data by gluing together compactly 
supported solutions: everything works because what we end up with is a locally finite sum. 
Step 2. Here, we describe in detail the localization argument. 

We introduce a disjoint covering |j£T Qi °f almost all of M'^, where {Qi} are open d- 
dimensional cubes with unit edge length. Let u as in fl2.ip . then u = ^^'Uj, where 
Ui = uxQi and XQi denotes the characteristic function of the set Qi. In Step 3 we construct 
a solution of the Cauchy problem 

dtUi{t, x) + div [b{t, x)ui{t, x)) =0 
Ui{0,x) = Ui{x) 

satisfying Ui G C°([0, +oo[; L\R'^) - s) and 

supp Ui C [0, T] X + Blt{0), (2.15) 

where L = \\b\\L°°. We set 

+ 00 

u{t,-) ■= ^^Uiit,-). 

1=1 

By noticing that the transport equation is linear and that the previous sum is locally finite 
by (I235D . we get that u G C°([0, +oo[; Ll^iM'^) - s) is a solution of flO) . 

Step 3. We fix T > and we construct a strongly continuous solution Ui of (12.141) defined 
for t G [0,T] and satisfying (12.151) . To simplify the notation, we write u, u and Q instead of 
Mj, Ui and Qi. 

We first approximate p and 6 by := p * rj^ and := (bp) * rj^/p^ respectively. Here rj^ is 
a sequence of smooth convolution kernels: the function p^ is always strictly positive by fl2.9l) 
and hence be is well defined. By using the notation L = we get that all vector fields 

be are bounded by L. Also, we can assume that the support of u*'r]s is strictly contained in 
Q + B^{0). Because of (12. 9p . there exists a constant M satisfying 

< Pe{t, y) < M for every (t, y) such that t G [0, T] and x G Q + fii+3LT(0). (2.16) 

Let denote the flow generated by b^ and let denote the inverse of Fg, i.e. we have 
Zeit, Ys{t, y)) = y and Ye{t, Z^it, x)) = x. Then 

\Ye{t, y)-y\< llfeelU-t < Lt for every [t, y) G [0, T] x 

and analogously 

\Z,{t,x)-x\< Lt for every (t, x) G [0, T] x W^. 

Let JYeit.x) denote the Jacobian of Fe(t,x): since ps(t,Ys{t,y))JYir{t,y) = ps{0,y), by 
relying on (I2.16P we deduce the bound 

^ < JYe{t,y) < for every t e [0,T], y e Q + B^+2Lt{0). (2.17) 
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Since X E Q + i?i+LT(0) implies Z^{t,x) E Q + -Bi+2lt(0) we obtain 

^ < JZ,{t,x) < M\ for every t G [0,r], x G Q + B^+lt{^). (2.18) 

We now define := u * rj^^/pi^iQ,-). Then, by setting le{t,x) := li.{Z^{t,x)), we get a 
solution of the transport equation 

dtle + 6, ■ V/, = 

Since ||Ae||z,oo < ||Ae||Loo < MUmH/^oo, {\} is weakly-* compact in L°°. Moreover, observing 
that peit, ■) pit, ■) strongly in L\^^{W^) for every t > 0, up to subsequences (not explicitly 
labelled for simplicity of notation) we get 

Xepe ^ u weakly-* in ( [0, T] x ) 

for some bounded function u whose support is contained in [0,T] x [Q + Blt{0))- First we 
point out that, since X^Pe solves 

dt{XePe) + div (beXepe) = , 

then u solves 

dtu + div(6n) = 

in the sense of distributions. Moreover, by extracting a further subsequence we can assume 
that 

XePeit, ■) ^ u{t, ■) weakly-* in L~ (R'^) for every t G [0, T] n Q . (2.19) 

Since Xi;pi;{t,-) is weakly-* continuous, uniformly in e (see for instance [71 Lemma 8.1.2]), 
fl2.19p ensures that 

XePeit, ■) ^ uit, •) weakly-* in L°° (M^) for every t E [0, T] . 

Thus, the only issue we are left to address is that uit, ■) converges strongly in Ll^^iM.^) to 
m(0, ■) = u when t — > 0^. By the lower semicontinuity of the norm with respect to the 
weak-* convergence in L'^ 



' oo 



|n(t, a;) — u(x)| < liminf / \Xepeit, x) — Xeix)peiO, x) \ dx. (2.20) 
Now, for any given e < 1, we exploit the decomposition 

\XePeit,x) - Xeix)peiO,x)\dx < / Peit, x)\leit, x) - Xeix)\ dx 

(2.21) 

+ \Xeix)\\peit,x) - peiO,x)\dx. 
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We first focus on tlie first term in (12.211) : for any G C^(]R'') we liave 



\le{t,x) - Ae(x)| dx 



< 



\Xe{Ze{t,x)) - Xe{x) \ dx 

Q+Bi+Lt{0) 

UZe{t,x))-tlj{Z,{t,x))\dx 



Q+Bi+Lt{0) 



[Z,{t,x))-^{x)\dx + 



Q+Bi+LtiO) 



(x) — Ae(x)| dx 



Q+Bi+LtiO) 



< 



Q+Bi+Lt(0) 



\\{Z,{t,x)) - ^{Z,{t,x))\ dx + Lip(^)^'^(g + 5i+it(0)) Lt 



By exploiting tlie bound (12.181) on JZ^ we tlien get 



(2.22) 



\K{Z,{t. X)) - i,{Z,{t,x))\dx < M'Wij - KWm^.y (2.23) 

'Q+Bi+Lt{0) 

Hence, by combining (I2.16p . (12.221) . (12.231) and taking into account tlie strong convergence 
of Ae to u/p{0, ■) as 5 ^ 0, we eventually obtain 

limsup / pi;{t, x)\l!;{t, x) — Xe{x)\ dx 



< M(l + 



u 



P(0, 



+ MU^{'ilj)^\Q + B^+lM) Lt. (2.24) 



We now focus on the second term in (I2.2ip : since A^ is bounded by M||m||lcx), by relying on 
the properties of the convolution we have 



|A£(x)| |pe(t, x) — Pe(0, x)\dx < M\\u\\L'y= I |p(t, x) — p(0, x)\ dx < M\\u\\LoaUo{t) . 

Jq+Bi{0) 

(2.25) 

In the previous expression, u denotes a modulus of continuity u : [0, +oo[^ [0, +oo[ for the 
map t I— > p(t, ■)xq+Bi(o)) which is strongly continuous in L-'^(M"'). 

Finally, we combine fl2:20D . ([MID, (12:251) and we eventually obtain 

|n(t, x) - u{x)\dx <M(1 + M')||^ - n/p(0, OIIl^m^) 

+ MLip(^/;)^^(Q + Bi+Lt{0)) Lt + M\\u\\L^u{t). 



By letting t I and using the arbitrariness of ip: we deduce that t u(t, ■) is strongly 
continuous at in L^. This concludes the proof of Theorem 12.71 □ 



WELL-POSEDNESS OF CONTINUITY EQUATIONS AND APPLICATIONS 



15 



We point out that in the proof of Theorem 12.71 the assumption that b and p are both 
BV\oc in space and time is only needed in the Uniqueness Part to show that distributional 
solutions are renormalized on ]0, +cx)[xR'^ and in the first step of the Existence Part to 
apply Lemma [2TT] (ii) and get that the map t i— >■ u{t, •) is strongly continuous in Ll^^(M.'^) at 
any r > 0. Hence, if we additionally assume that the divergence of 6 is a locally integrable 
function, we can relax the regularity assumption on both b and p. More precisely, assume 
that b E Ll^^(jO, +oo[; BVioci'^'^;'^'^)) is a bounded, locally strongly nearly incompressible 
vector field such that div 6 G Lj'j^^(]0, +cx)[x]R'^). Then we can directly apply Ambrosio's 
theorem [1] to deduce that distributional solutions are renormalized. By arguing as in the 
proof of Lemma l2.1l (ii). we obtain that, if u is a locally bounded solution of (12. ip . then 
t I— > u{t, ■) is strongly continuous from ]0, +oo[ in Lj'Qj,(M'^). By combining all the previous 
considerations, we get: 

Theorem 2.8. Let b E Ll^^(jO, +oo[; BViocO^'^',^'^)) bounded, locally strongly nearly 

incompressible vector field such that div 6 E Ll^^(JO, +oo[xM.'^Y Then there exists a locally 
bounded solution u E (^[0, +oo[; Ll^^(M'^) — s) to the Cauchy problem for the continuity 
equation (12. ip . Furthermore the solution is unique in this class. 



3. An example of nonuniqueness 

Roughly speaking, the results presented in Section [231 allow to weaken the assumptions on 
the summability of the BV norm of the vector field, at the price of a restriction of the class 
of solutions considered. In the model case of a divergence free vector field (which is trivially 
locally strongly nearly incompressible), comparing Theorem 12.81 and Ambrosio's theorem [T] 
we see that we swap summability of the BV norm up to t = for the condition of strong 
continuity of the solution. 



3.1. Depauw's vector field. In [21] Depauw constructs an example of nonuniqueness for 
the transport equation relative to a bounded time-dependent divergence free vector field in 
the plane. More precisely, Depauw exhibits a vector field a(t,x) E L°°([0, 1] x R^;M^), with 

diva = and a G LJ^^QO, 1]; 5Moc(M^; M^)) , 
such that the Cauchy problem 



dtU + div (an) 
m(0,x) = 



(3.1) 



has a nontrivial solution n ^ 0. Observe that, being divergence free, the vector field a is 
locally strongly nearly incompressible (see Definition 12.60 : simply take p{t,x) = 1. Thus, 
Theorem 12.71 implies that u cannot be strongly continuous in time, since the constant 
trivially provides the unique strongly continuous solution of (13. ID . 
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We now want to briefly sketch the construction of [21]. Let us first consider a dyadic 
subdivision (up to neghgible sets) of the time interval, i.e. 

oo 

[0, 1] = U 4 with 4 = 

k=l 

For t & Ik the vector field a{t, ■) does not vary in time, and consists of "vortexes" on a 
pattern with scale 1/2^, arranged in such a way that the divergence is zero and, for any 
compact set K cM"^, 

Ht,-)\\BViK)-'^' fortG/fc. (3.2) 
This allows to construct a solution which undergoes, on each interval 4, a more and more 
refined mixing, as t — > 0"^. More precisely, going backward in time on each interval Ik, the 
solution is rearranged from a function which oscillates between ±1 on a chessboard of size 
1/2^^^^ into a function which oscillates between ±1 on a chessboard of size 1/2^^. This in 
particular implies that 

u{t,-)^0 in L°°(R2) - w*, 

but 

which is coherent with Theorem 12.71 

3.2. Approximability of the solution with smooth maps. The proof of Theorem 12.71 
directly shows that any limit of solutions with smooth approximating vector fields IS, as a 
matter of fact, strongly continuous. Hence we get that the solution u constructed in [21] 
cannot be constructed by approximation. 

It is also worth to mention a connection with the results in [10], regarding the density 
of smooth functions in the space of the solutions of the transport equation. Namely, Theo- 
rem 2.1 in [To] asserts that the approximability of a solution with smooth functions is equiv- 
alent to the uniqueness property for both the forward and the backward Cauchy problems. 
We note that, while the vector field of [2T] provides a counterexample to forward uniqueness, 
the same vector field enjoys uniqueness for the backward Cauchy problem. More comments 
and variations on the construction by Depauw are also presented in jlOj . 

3.3. Strong continuity of the vector field does not imply uniqueness. The vector 
field described in Section 13.11 provides a counterexample to the uniqueness in the class of 
weakly continuous solutions, but 

a^C°([0,l];LL(M2;M2)). 

Thus, we could wonder whether Theorem 12.71 holds in the larger class of weakly continuous 
solutions, provided we impose strong continuity of the vector field, i.e. 

beC%[0,+oc[;LU<^'))- 
This means that we are led to the following question. 



1 1 
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Question 3.1. Let b G -BVioc(]0, +oo[xR''; R'')) be a bounded, locally strongly nearly in- 
compressible vector field. Assume that an admissible density p in the definition of locally 
strongly nearly incompressibility belongs to i?Vioc(]0, +oo[xR"'). Assume moreover that 

Is it true that the Cauchy problem fl2.ll) has a unique locally bounded solution u in the class 
CO([0,+oo[;L~(M'^)-w;*)? 

Notice that, if the answer to Question 13.11 were positive, we would be able to prove one 
of the results regarding the well-posedness of the chromatography system, the one presented 
in Theorem 14.81 in the more general class of solutions 

U = {uuU2) G C%[0, +oo[; LZi^; R^) - w*) . 

However, we can show the following negative result: 

Proposition 3.2. The answer to Question \ 3.1\ is negative. More precisely, there exists a 
vector field c(t,x) G L°°([0,T] x R^;]R^) satisfying the following properties: 

(i) CG Wioc(]0,T] XR2;R2); 

(ii) dive = (and thus c is in particular locally strongly nearly incompressible); 

(iii) ceC^{[0,T];L°°{R^;R^) - s); 

(iv) The Cauchy problem 

dtV + div (cv) = 

v{0,x) = ^ ' ' 

has a bounded solution v ^ 0, with v ^ ^([0, T]; Lj'Q^(R^) — s) . 

Proof. We briefiy illustrate how it is possible to modify the construction in [21] in such a way 
that the resulting vector field is strongly continuous. The vector field a is the one described 
in Section [3TT1 

Step 1. Strong continuity at time t = 0. The vector field a has length between 
and 1 in all the regions in which it is nonzero. This length is chosen in such a way that the 
solution is mixed from a chessboard of size 1/2^^^ into a chessboard of size 1/2^ precisely in 
time 1/2''. 

To construct the vector field c, we decrease the strength of the vector field in each interval 
in such a way that the discontinuity at t = is ruled out. First, we substitute each interval 
Ik with an interval Jk such that 

We notice that Yl^=i-^^iJk) = T < +oo, thus we are again concerned with a finite interval 
of time. In each interval we then simply multiply the vector field a relative to the interval 
Ik by 1/k. Thus, the solution undergoes the same mixing, but in a time which is k times 
larger than the one in Depauw's example. Notice that we precisely have ^^(J^) = k^^{Ik)- 
In this way, we have the convergence c(t, ■) ^ strongly in L°°(R^) as t | 0. 
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Step 2. Blow up of the BV norm. Let us check that the spatial BV norm of the 
vector field c{t, ■) blows up as t J, 0. We first notice that for any compact set K C we 
have 

||a(t,-)||w{AO ~ 2*^ fort eh, 
thus \\a{t, ■)\\BViK) ~ In our case, we simply have 

1 

Noticing that 

' — • . oo ^ 

2J - ^ 2J 



\c{tr)\\ 



BV{K) 



— for t G Jh 



E 

j=k+i 



k_ 
2fc 



we have that t > when t e Jk, and this implies 



Ht,-)\\ 



BV(K) 



>1 
~ t 



^)), which is coherent with the 



Let us notice that in each 
However, at the extrema of 



This shows that c does not belong to L-*^ ([0, T]; i?Vioc(IR" 
uniqueness result of [I]. 

Step 3. Strong continuity at any time t e [0,T] 
interval Jk the vector field c is constant with respect to time 
the intervals Jk the vector field c has a jump discontinuity. To rule out this phenomenon, 
we choose for every time interval an "activation function" tpk{t), and we substitute the 
vector field c(t, x) previously described with c{t, x)ijjkit), for t G Jk- We only need to choose, 
on each Jk, the function ipk satisfying 



< 2 and 



Ms)ds = ^\Jk 



Jk 



In this way, the solution undergoes the same mixing as before in each time interval Jk, but 
the discontinuities at the extrema of the intervals Jk are ruled out. It is also immediate that 
dtc is a locally finite measure, hence c G BVioc{]0,T] x M*^). 

Step 4. Conclusion. In this way, we have constructed a vector field c(t, x) e L°°{[0, T] x 
M^;M^) satisfying properties (i)-(ii)-(iii). The nontrivial solution v of (13.31) has the same be- 
haviour of the nontrivial solution u of (13.11) described in Section 13.11 at the Lagrangian level 
we have simply reparameterized (with the same function) all trajectories, "stretching" the 
original vector field a in such a way that we gain strong continuity of the vector field c with 
respect to the time. □ 



4. Applications to the chromatography system 

4.1. The chromatography system and classification of the entropies. In this section 
we are concerned with the system (11.11) coupled with the initial conditions 



Mi(0, x) 



Ui{X) 



U2{0,x) 



U2{Xj 
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This system belongs to the Temple class [33] and arises in the study of two component 
chromatography. Here Ui, U2 : [0,+oo[x]R M are nonnegative functions which express 
transformations of the concentrations of two solutes, see [121 page 102]. Here we make some 
preliminary considerations by introducing a change of variables and discussing a related clas- 
sification of the entropies. Our well-posedness results for (II. ip are discussed in Sections 14.21 
and I 



By introducing the change of variables 

V = Ui+U2, 



W 



U1-U2, 



(4.2) 



system (II. ID becomes 



dt'u + 9^ 



dtw + 



(4.3) 



w 



0. 



Namely, (II. ip splits in the coupling between a scalar one-dimensional conservation law and a 
linear continuity equation. Note that the coefficient of the second equation in (14. 3 p depends 
on the solution of the first equation, but the first equation in (14. 3 p does not depend on the 
solution of the second one. 

This decoupling allows for a classification of all the entropies of system (11.10 : 

Lemma 4.1. Let (77, g) : ^ M x M be twice continuously differentiable. Then {ri,q) is an 
entropy- entropy flux pair for system (II. ID (in the sense of Definition \A . i\) if and only if the 
following holds: there exists an entropy- entropy flux pair (f/, g) : M ^ M x M for 

V 



dtv + 
satisfying 

r]{ui, U2) = fj{ui + U2) + C{ui - U2) , 







q[Ui,U2) = q{ui + U2) + 



C{ui - U2) 



l + Ui 

for some real constant C . Moreover, rj is convex if and only if fj is convex. 



U2 



(4.4) 



(4.5) 



Proof. First, observe that (77, q) is an entropy-entropy fiux pair for system ( 11. ID if and only 
if 

fv + W V -w^ 
{Ti,q){v,w) = {ri,q) ( 

is an entropy-entropy fiux pair for system (14. 3p . Thus, in the following we focus on system 
(14. 3p . having fiux given by 



G{v,w) 



1 +v 
w 

\TT~v J 
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Being t] an entropy, we have the compatibihty condition 

= 0. 



Since 



curl ( V r]{v, w) ■ DG{v, w] 
/ 1 



(4.6) 



DG{v,w) 



\ 



'1 



-w 



\ (1 



we get from (14. 6 p 



namely 



curl 



(1 ' l + v 

vd^yj-q + wdyj^T] = , 



By setting 7 = d^rj, we obtain 

V7(t', w) ■ {v, w) = for any vector {v, w). 

Hence, 7 is constant along any ray departing from the origin, and since 7 is of class C^, it 
has to be constant on the whole M^. Thus, 

ri{v, w) = fj{v) + Cw 

for some real constant C. Being in the scalar case, fj is trivially an entropy for (14. 4p : the 
flux is defined by 

q\v) = fi\v)g\v) , 

where g{v) = v/{l +v) is the flux of (14. 4p . Moreover 77 is convex in (w, w) if and only if fj is 
convex in v. Hence, up to constants, the entropy flux q satisfies 

Cw 



q[v,w) 



q{v) 



1 



This concludes the proof of the first implication: any entropy-entropy flux pair [rj, q) for 
f lLT]) satisfies fli3|) . 

Conversely, a straightforward computation ensures that, for any entropy-entropy flux pair 
{fj,q) for (14. 4p . (14.50 deflnes an entropy-entropy flux pair (r/, g) for (II. ip . □ 

Heuristically, Lemma 14.11 can be interpreted as follows. By introducing the change of 
variables (14. 2 p we select a direction w in the phase space (wi, M2), along which the solution is 
simply transported. Consequently, the entropy dissipation in the direction w is zero. Hence, 
imposing that {ui,U2) is an entropy admissible solution of (II. ip means imposing that the 
entropy dissipates along the orthogonal direction v. This interpretation is coherent with the 
ideas underlying the analysis in [9]. 

By relying on Lemma 14. we obtain a one-to-one correspondence between entropy ad- 
missible solutions of the system (ILip and distributional solutions of system (14. 3p such that 
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the first component v is an entropy admissible solution of tlie scalar conservation law fl4.4p . 
Indeed the following holds: 

Proposition 4.2. Let ui, U2 G Li'^^(R). Then (1x1,^2) G Li'^j,([0, +oo[xM; M^) is an entropy 
admissible solution of (11. ip satisfying the initial condition 

(mi(0,-),W2(0,-)) =(wi,W2) (4.7) 
if and only if [v, w) = {ui + U2, Ui — U2) verifies the following conditions: 
(1) V is an entropy admissible solution of the Cauchy problem 



dtv + dj^] =0 



(4. 



.1^(0, ■) = ui + U2; 
(2) w is a solution in the sense of distributions of 

dtw + d^' ^ 



1 



(4.9) 



w{0, ■) = ui- U2. 

The initial data in (14.81) and (14. 9 p are assumed in the same sense (L\^^ — s or — w*) as 
the datum (^1,^2) in (14. 7p . 

Proof. By definition {ui,U2) is an entropy admissible solution of (11.11) if and only if 

dtviui, U2) + d^q{ui, U2) < (4.10) 

in the sense of distributions for any entropy-entropy flux (77, q) with rj convex. By applying 
f l45ll . (141011 becomes 

C{ui - U2) 



fj{ui + U2) + C{ui - U2) 



q[ui + U2) 



<0. (4.11) 



I +Ui + U2_ 

for any entropy-entropy flux pair {fj,q), with fj convex. Since (14. lOp is assumed to hold for 
any entropy-entropy flux (77, q) with rj convex, in (14. lip we can take any C G M. Thus, 
requiring (14. lip is equivalent to imposing that w = Ui — U2 is a. solution in the sense of 
distributions of (14. 9 p and v = ui + U2 is an entropy admissible solution of (14.81) . 

Finally, as the change of variables (14. 2 p is linear, the initial data in (14.81) and (14.91) are 
assumed in the same sense as the datum {ui,U2) in (14. 7p . □ 

4.2. Well-posedness results in the class of weakly continuous functions. We first 
apply Theorem 12. 5^ obtaining a well-posedness result under the following assumptions on 
the initial datum U = {ui, ^2): both ui and U2 are nonnegative (which is the physical range), 
locally bounded, and the sum U1 + U2 has locally bounded (but possibly large) total variation. 
We introduce the following set: 
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Definition 4.3. We denote by T the set of functions U = {ui, U2) : M ^ M x M such that: 

(.Fl) Ui, U2 e (R); 

(^2) ui > 0, U2 > a.e. in M; 

We are now ready to introduce our result: 

Theorem 4.4. Let U = {ui,U2) G JF. Then there exists a unique solution 

U = (Ml, M2) e C° ([0, +oo[; L~ (M; R') - w*) 

in the sense of distributions of the Cauchy problem (Il.ip - fl4.1l) . among those with values in 
T and satisfying the entropy admissibility condition given by Definition \A.l[ 
In addition, this unique solution U satisfies the regularity property 

f/GC°([0,+oo[;LL(M;M2)_s). 

Proof. By applying the change of variables fl4.2l) and by relying on Proposition 14.21 we reduce 
the problem to the study of the Cauchy problems (14.81) and (14. 9p . We then proceed in several 
steps. 

Step 1 . By applying Theorem IA.2I and exploiting Remark I A. 71 we obtain that the first 
equation in (14. 8 p admits a locally bounded entropy admissible solution v G C*^ ([0, +oo[; Lj^^(R) 
w*y Moreover, thanks to estimate (lA.8p . for every t > we have v(t,x) > for a.e. x G M. 
Estimate (]A.9P ensures that v G L^^[[0, +oo[; i?V[oc(M)). Finally, by relying on (14.80 we get 
that V has also bounded variation with respect to the time, thus v G -BVioc([0, +oo[xR). 

Step 2. We now apply Theorem 12.51 with 

b(t,x) = 7 r 

and p = V, deducing existence and uniqueness for the solution w of ( 14.90 in the class 

C'{[0,+CK>[-L^M)-w*)nL^{v), 

where L°°{v) is defined as in Definition 12.31 Furthermore Theorem I2.5( iii) implies that 
\w(t,x)\ < v{t,x) for a.e. x G M, for all t G [0,+oo[. 

To show that the solutions ui and U2 are nonnegative, it suffices to observe that we have 
ui = [v + w)/2 and U2 = {v — w)/2, and we have just shown that \w\ < v. This completes 
the existence part. 

Step 3. We want to show that the solutions v{t, ■) and w{t, ■) constructed above are 
strongly continuous with values in Ll^^(M.'^). The strong continuity of t 1— >■ f (t, ■) is a con- 
sequence of Theorem IA.5I (since v i— > v/{l + v) is uniformly concave we can exploit the 
considerations in Remark lA. 60 . Concerning w, it suffices to apply again Lemma [2.1( 11) with 
q = V (which has just been shown to be strongly continuous) and u = w/v (see Remark l2.4p . 

Step 4. We finally show uniqueness. By the uniform concavity of the flux function 
V v/{l + v) in (14.80 . Theorem I A . 51 implies that every weakly continuous entropy admissible 
solution of (14. 8 p is in fact strongly continuous at t = 0. Hence the uniqueness of v follows 
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from Theorem IA.2[ Concerning i/;, it is sufficient to observe that, since condition {J^2) gives 
Ui and U2 are nonnegative, we have 

\w\ = |mi — M2I < Ui + U2 = V . 

Thus w G L°°{v), and we get uniqueness from Step 2. □ 

As an immediate consequence of Theorem 14.41 we obtain: 

Corollary 4.5. There exists a unique semigroup 

St : [0,+cx)[ X J^-> 

satisfying: 

(i) St is continuous with respect to t in Lj^j,(M;M^) endowed with the weak-* topology; 

(ii) For every If = {ui,U2) G ^, StU provides a distributional solution of the Cauchy 
problem ( ll.ip - (]4.ip satisfying the entropy admissibility condition given by Defini- 
tion \A.li 

Remark 4.6. By relying on the proof of Theorem 14.41 one actually gets a slightly sharper 
result than the one given in the statement of the theorem. Namely, for any U = {ui, U2) G 
there exists a unique 

f/ = (ni, U2) e C%[0, +oo[; L~ (M; R') - w*) 

solution in the sense of distributions of the Cauchy problem fll.ll) - fl4.ip . satisfying the entropy 
admissibility condition given by Definition lA.ll and such that U{t,-) = ■), ■U2(^, ■)) 

verifies conditions (J-'l) and {^2) in Definition 14. 3[ For such a solution, property is 
automatically satisfied by U{t, ■) = [ui(t, ■), U2(t, ■)) ^oi every t > 0. 

4.3. Well-posedness results in the class of strongly continuous functions. We now 

apply Theorem 12. 71 obtaining well-posedness results for ( ]l.ip -( l4rTI) under assumptions on the 
initial data {ui,U2) different from those considered in Section I4l2l and in a different class of 
solutions. Namely, we relax the regularity assumptions imposing only that [ui, U2) G LJ^j,(M). 
The price we have to pay is requiring that the sum ui + U2 is well-separated from on any 
compact set, in the sense of fl4.12p . 
We ffist introduce the following set: 

Definition 4.7. We denote by Q the set of functions U = (mi,M2) : M — R x R such that 

(gi) Ml, U2 G (R); 

(^2) Ml > 0, M2 > a.e. in M; 

(^3) for every R > there exists Sr > such that 

ui + U2>6r a.e. in] - R,R[. (4.12) 

We are now ready to state our result: 
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Theorem 4.8. Let U = {ui,U2) G Q. Then there exists a unique 

solution in the sense of distributions of the Cauchy problem (Il.ll) - fl4.ip and satisfying the 
entropy admissibility condition given by Definition \A.l\ In addition, the unique solution U 
takes its values in Q. 

Proof. As in the proof of Theorem 14.41 we reduce to study the Cauchy problems fl4.8l) 
and (14.91) . We then proceed in several steps. 

Step 1. By applying Theorem I A. 2 1 and exploiting Remark I A. 71 (with f{z) = z/{l + z)) 
we obtain that (14. 8 p admits a nonnegative entropy admissible solution v G Lj^^([0, +cxd[xM). 
Observing that constant functions are entropy admissible solutions of (14. Sp . thanks to esti- 
mate (lA.SP and (14.120 we easily obtain that, for every R, t > 0, 

v{t, x) > Sn+Lt for a.e. x e] - R, R[. (4.13) 

Proposition IA.4I combined with Remark IA.6I ensures that 

^GL~(]0,+oo[;Wioc(M)), 

and using (14. 8 p we deduce that 

V G 5Moc(]0,+oo[xM) . 

The fact that 

^;GC°([0,+oo[;LL(M)-s) (4.14) 

is a consequence of Theorem IA.51 thanks to the uniform concavity oi v ^ v/{l + v). 
Step 2. Let 

h{t,x)= ^ (4.15) 
1 + f (r, x) 

where v is the same function as in Step 1. Then 6 is a bounded vector field belonging to 
-BVioc(]0, +cxd[x]R). Moreover b is locally strongly nearly incompressible: by taking p = v 
in ( 12. 9p and (I2.10p we get that Definition 12.61 is satisfied, as v solves 

dtv + d,{bv) = 

and verifies (I4.13P and (14.141) . Then, Theorem 12.71 ensures that the Cauchy problem (14.90 
admits a unique solution w G C°([0, +oo[; LlM) - ^^)- 

Step 3. Let b be the same vector field as in (I4.15p . Then v and w are the unique solutions 
(by Theorem 12. 7p of 

d,z + d {bz) =0 ^ ^ _ ^^^^g^ 

2;(U, x) = z[x) , 

with z = V and z = w respectively. We claim that the Cauchy problem (I4.16P satisfies the 
following comparison principle: 

z{x) > for a.e. x G M =^ for every t > 0, z{t,x) > for a.e. x G M. (4.17) 
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This implication is proved in Step 4. To conclude with the existence part, we apply (14.171) 
with z = V + w and z = v — w, noticing that at t = we have v + w = 2ui and v — w = 2u2, 
which are both positive by condition {Q2) in Definition 14.71 

Step 4. We now show implication (14.171) . As it comes from the existence part in the 
proof of Theorem 12.71 the unique solution of (I4.16P is obtained by approximating the vector 
field b with smooth vector fields be ^ b strongly in Lj^^.^^, +oo[x]R) and considering the 
solutions Ze of the corresponding problems 

dtZs + d^{beZe) = 

Ze{0,x) = z* r]e{x) . 

(With the notation of the proof of Theorem 12. 7[ z^ = lePe-) Since ^ > a.e. in M, for 
every t > we have Ze{t,x) > for a.e. x G M. Hence, observing that z^ converges to z 
weakly-* in lv[^^([0, +oo[xM), we conclude that z{t,x) > for a.e. {t,x) G [0, +oo[x]R. Since 
z G C°([0, +oo[; LlM) - imphcation (liTfll follows. 

Step 5. Uniqueness is a simple issue: since v is by assumption strongly continuous, it is 
unique by Theorem lA. 21 while w is unique thanks to the uniqueness part of Theorem 12. 71 □ 

Theorem 14.81 guarantees, in particular, that the domain Q is invariant for admissible solu- 
tions of (II. ip : if the initial datum (^1,^2) belongs to Q, then [ui(t, ■),U2(t, ■)) G Q for every 
t > 0. The following is an immediate consequence of Theorem 14. 4t 

Corollary 4.9. There exists a unique semigroup 

St : [o,+cx)[ xg^g 

satisfying: 

(i) St is continuous with respect to t in L\^^{M.]'E?) endowed with the strong topology; 

(ii) For every U = {ui,U2) G g, StU provides a locally hounded distributional solution of 
the Cauchy problem (Il.ip - (l4.ip satisfying the entropy admissibility condition given by 
Definition \A.1\ 

Remark 4.10 (Extension to the k x k system). The k x k chromatography system is the 
following: 

dtUi + d,j,(— — ■ )=0, i = l,...,fc. 

\1 + Ui + . . . + Uk J 

By doing the change of variable v = Ui + . . . + Uk, Wi = Mj+i for i = 1, . . . , k — 1, our system 
becomes 

dtv + dj^^=0 

dtWi + dr, { ] = , i = l,...,k-l. 



1 + 

All the results proved in this section, including the classification of the entropies-entropy 
fiuxes given in Lemma [4.11 readily extend to this case. 
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5. Applications to the Keyfitz and Kranzer system 

As a byproduct of our analysis, in this section we apply the results discussed in Section [2] 
to the Keyfitz and Kranzer system 



^ ^ (5.1) 



^ U{0,x) = U{x). 

The function U takes values in and depends on {t, x) G [0, +oo[ xR''. For each a = 1, . . . ,d 
the function /q, : R ^ M is smooth. As mentioned in the introduction, to study (15. ip one 
can first obtain the modulus of the solution p = \U\ by solving 



dtp + div{f{p)p)=0 
p(0,x) = |f/|(x). 



(5.2) 



where / = (/i,...,/^). Then, one gets the solution U = {p9i, . . . , pOk) by solving the 
continuity equations for the components of the "angular" part 

dt{pei)+dw{fip)p9i)=0, t = l,...,k. (5.3) 

This strategy has been exploited by Ambrosio and De Lellis [5] and by Ambrosio, Bouchut 
and De Lellis [3], and was inspired by considerations in Bressan [13]. Before discussing 
what we get by combining this strategy with Theorems 12.51 and 12.71 we need to provide the 
following definition: 

Definition 5.1. Let U he a. locally bounded function solving (15. ip in the sense of distribu- 
tions. Then U is a renormalized entropy solution if p = \U\ is an entropy admissible solution 
of (15. 2p (in the sense of Definition lA.ip such that limt^o+ p(t^ ') = 1^1 the strong topology 



of 



As pointed out in Ambrosio, Bouchut and De Lellis [3] by relying on a classification of 
the entropies for (15.10 due to Frid [21], under quite general assumptions on the function / 
any renormalized entropy solution is indeed an entropy admissible solution of (15. ip (see for 
example Dafermos [HI Chapter IV] for the definition of entropy admissible solution of a 
system of conservation laws in several space dimensions: this notion is the multidimensional 
analogue of Definition lA.ip . A complete proof of this implication can be also found in the 
notes by De Lellis [201 Proposition 5.7]. Conversely, an example discussed in Bressan [131 
Section 3] shows that, even in the one-dimensional case, in general there might be entropy 
admissible solutions of (15. ip that are not renormalized entropy solutions. The same example 
shows that, in general, entropy admissible solutions of (15. ip are not unique. 

Existence, uniqueness and stability results for renormalized entropy solutions of (15. ip 
were obtained in Ambrosio and De Lellis [5] and Ambrosio, Bouchut and De Lellis [3]. In 
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particular, in [3] the Cauchy datum U in (15. ip satisfies |f/| G L°° fl BVioc and \U\ can attain 
the value 0. 

By applying Theorem 12.71 we manage to relax the assumption \U\ G L°° fl -BVioc by 
requiring only |f/| G L°°. However, the price we have to pay is that we restrict to the one- 
dimensional case, we assume that the map p ^ f{p)p is uniformly convex, and we impose 
that the initial datum \U\ is well separated from 0. Similar results were already known: for 
example in Freistiihler [23] and Panov [30] well-posedness theorems in one space dimension 
for L°° initial data were obtained under more general assumptions than those we consider 
here. However, Theorem 15.21 below quickly follows from Theorem 12.71 so for completeness 
we provide the details of the proof. 

Theorem 5.2. Assume that the following conditions are satisfied: 

(1) there exists c > such that [f{p)p]" > c on M; 

(2) U G L°°(M;M^'); 

(3) For every R > 0, there exists Cr> such that 

\U{x)\>^>0 for a.e. X e[-R,R]. 

Then the Cauchy problem 

dtU + d,[f{\U\)U] =0 

(5.4) 

f/(0,x) = U{x) 
admits a unique renormalized entropy solution 

U G C°([0,+oo[; Lf„,(M,R^) - s) n ([0, +cx)[xM; R^) . 

Proof. Existence. The proof exploits the splitting of (15.11) in the coupling between (15.21) 
and (15.31) . We proceed in several steps. 

Step 1. We first construct the modulus p of the solution. Indeed, Theorem IA.2I ensures 
that there exists a unique locally bounded entropy admissible solution p to 

dtp + dt{fip)p)=0 (5.5) 

which satisfies 

hm p(t, ■) = strongly in L^^e W- (5-6) 

Step 2. Our goal is now applying Theorem 12.71 Let us check that all the hypotheses are 
satisfied. We first show that 

bit,x)= [f{p)]{t,x) (5.7) 

is a locally strongly nearly incompressible vector field, in the sense of Definition 12. 6[ By 
combining (lA.SP and conditions (2) and (3) in the statement of the theorem we get that for 
every R> and T > there exists a constant Cr^t > such that 

^ < p{t,x) < \\U\\loo for a.e. (t,x) G [0,T] x 5^(0). (5.8) 



C 



R,T 
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Also, by Theorem IA.5I and by assumption (1), the map t t— > p{t, ■) is strongly continuous 
with values in Lj'Q^(M). 

Thanks to Proposition lA.4l and assumption (1), we deduce that p G Lj^^(]0, +oo[; -BVioc(K))• 
Moreover, by exploiting the equation satisfied by p we also have p G BViocQO, +oo[xM.), and 
since / is smooth, the same regularity is inherited by the vector field b defined in ( 15. 7p . 

We can then apply Theorem 12 . 71 and conclude that, for every U = {ui, . . . , Uk) G L°°(]R; M^') 
the Cauchy problems 

dtUi + dx{bui) = 



u,{0,x)=u,{x), (5.9) 

admit a unique locally bounded solution U = {ui, . . . , u^) G C°([0, +oo[; ^^^^^^(IR; M'^) — s). 

Step 3. To conclude the proof of the existence part, we need to show that, for every 
t > 0, 

p{t,x) = \U\{t,x) fora.e. xgM. (5.10) 
Thanks to (15.81) . for every i = 1, . . . , k we can define 6i := Ui/p, and we get 

dt{pei)+d^{bpe,) =0, 

in the sense of distributions on ]0,+oo[x]R. Since the functions 6i are locally bounded, we 
can apply Lemma [2.1( i) on the time interval [r, +oo[ for an arbitrary r > to deduce that 

dt{pe^)+d,{bpe^)=o, 

in the sense of distributions on ]0, +oo[xR, or equivalently 



^* Vj^^' W ' =0. (5.11) 



Summing over i = 1, . . . , k the equations in ( 15. lip , we obtain 



P / VP 
Moreover, (15.90 and (15.60 imply that 

,5.13) 

p(o,-) 

Recalhng (15.70 and comparing (15.50 and ( 15.120 . we see that p and \U\'^/p are strongly con- 
tinuous in time, and solve the same Cauchy problem (the initial data coincide because of 
( I5.13P ). Thus, by applying Theorem 12.71 we obtain that, for every t > 0, p(t,x) = \U\(t,x) 
for a.e. x G M, as desired. 

Uniqueness. Let Ua and Ub be two bounded strongly continuous renormalized entropy 
solutions of (15. 4p . Then pA = \Ua\ and ps = \Ub\ are two entropy admissible solutions 
of (15.20 for which limt^o+ Pa(^, ■) = ^^^t-*o+ Psit, ■) = \U\ strongly in Lj'Q^(M). By the 
uniqueness part in Theorem IA.21 we deduce that, for every t G [0, +oo[, 

PA(t,x) = pB{t,x) for a.e. x G M. 
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Let p denote the common value of pa and pB- We have 

dtUA,i + d., {f{p)UA,) = , dtUs, + {f{p)UB,) = for every t = l,...,k 

and 

UA,i{0, ■) = UB,i{0, ■) a.e. in M, for every i = 1, . . . , k. 
Theorem 12.71 thus gives that for every i = 1, . . . , k and all t G [0, +oo[ 

UAA{t,x) = UB,i(t,x) for a.e. x G M, 
as desired. □ 

Remark 5.3. In Theorem 15.21 we restrict to the one-dimensional case d = 1 because in the 
proof we need to apply Olemik's estimate flA.lOp to gain a regularizing effect (this is also 
the reason why we impose that the function p i-^ /(p)p is uniformly convex). 

Remark 5.4. By exploiting the same arguments as in the proof of Theorem 15.21 but applying 
Theorem 12.51 instead of Theorem 12.71 one obtains the existence and uniqueness result proven 
by Ambrosio, Bouchut and De Lellis in [3]. Namely, assume that \U\ G L°° fl BViod^'^)- 
Then (15.11) admits a unique renormalized entropy solution U G C°([0, +oo[; L°°{R'^; M*^) — w*) . 



Appendix A. Systems of conservation laws in one space dimension 

For completeness, in this short Appendix we go over some results that are used in the 
paper. Consider a system of conservation laws in one space dimension: 

dtU + d^[F{U)] =0, where (t,x) G [0,+oo[xM and f/ G M^ (A.l) 

with F : M'^ ^ M'^ smooth. For a general introduction to the subject, we refer for example 
to the books by Bressan [12], by Dafermos [19] and by Serre [32]. In the following, we focus 
on the Cauchy problem, assigning the initial condition 

U{0,x) = U{x). (A.2) 

It is known that, even if the initial datum U is very regular, in general there is no classical 
solution to (lA.ll) -( lAl2l) defined on the whole time interval t G [0, +oo[. Examples of solutions 



starting from a datum U G C°° and developing discontinuities in finite time are available even 
in the scalar case k = 1. It is thus natural to interpret (lA.ip in the sense of distributions, 
by requiring that t/ is a locally bounded measurable function satisfying 

/ / Udt<^ + F{U)d^<l> dxdt = for any $ G C,^(]0, +oo[ xM;M^). (A.3) 
Jo Jr'- 



However, in general a weak solution of the Cauchy problem ( ]A.ip -( lAl2l) is not unique. To 



select a unique solution, various admissibility conditions have been introduced, often moti- 
vated by physical considerations (see again Dafermos [I9]). Here we focus on the entropy 
admissibility condition. 
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Definition A.l. Let 77 : M*"' ^ M and g : M'^ ^ M be twice continuously differentiable 
functions. Tlie couple {r],q) is an entropy- entropy flux pair foi flA.ip if 

Vr){U)DF{U) = Vq{U) for every U G M^ (A.4) 

where DF denotes the Jacobian matrix of F. Let U : [0, +00 [xM — M'^ be a bounded, 
measurable function. Then U is an entropy admissible solution of (lA.ll) if the following 
holds: for any entropy-entropy flux pair (77, q) with i] convex, and for any (f) E (]0, +00 [xM) 
nonnegative, we have 

/•+00 /• 

r]{U) dt(j) + q{U) d^(j) dxdt > 0. (A.5) 

Let us now focus on a single conservation law in one space dimension 

dtu + dj{u) = 0. (A.6) 

Kruzkov [27] proved existence and uniqueness results for the Cauchy problem in the class of 
entropy admissible solutions (in the sense of Definition lA.ll) assuming the initial condition 
in strong sense: 



Theorem A. 2 (Kruzkov, [27]). Let f be a smooth function and assume that u E 
Then there exists a unique function u G L°°([0, +oo[xM) satisfying the following properties: 

(1) u is a solution in the sense of distributions of ( 1A.6[) ." 

(2) the solution u is entropy admissible; 

(3) lim uit, ■) = u in the L\ (M) topology. 

Furthermore, the solution u enjoys the following properties: 

Theorem A. 3 (Kruzkov, [27]). Let f be a smooth function, and for u, v E L°°(M) let u, 
V G ([0, +oo[xM) be the solutions provided by Theorem \A.^ to the Cauchy problems for 
( ]A.6|) with initial data u and v, respectively. Define further 

L := sup ||/'(z)| : |^| < max{||u||oo, ||^^||oo}|- (A.7) 
Then, for any R > and for all t > we have 

R j-R+Lt 

\u{t,x) — v{t,x)\ dx < [u{x) — v{x)~\ dx , (A.8) 

R J-R-Lt 

where [■]"'" denotes the positive part. Moreover, if we assume u G i?yioc(IR) and denote by 
ToiYsi R{u(t, ■)} the total variation of the function u{t, ■) on the interval ] — R, R[c M, then 
for all t G [0, +00 [ 

TotVar^{M(t, ■)} < ToWaTR+uiu} . (A.9) 

Under an assumption of uniform convexity on the flux, namely f"{u) > c > for any 
M G M, the classical Oleinik's result [29] asserts that the solution of the Cauchy problem 
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provided by Theorem IA.2[ with initial datum u G L°°(]R), enjoys the following one-sides 
estimate: for alH > it holds 

u{t,y) - u{t,x) <- — - for all X < y, X, y eR\Nt, (A.IO) 

with Nf Lebesgue negligible (possibly depending on t). In particular, flA.10|) implies that the 
space distributional derivative dxu{t, ■) of the solution u satisfies the upper bound 

dxu{t, ■) < — in V'{R), for all t > 0, (A.ll) 

ct 

and this implies in particular that dxu{t, ■) is a measure. The following fact is well known, 
but we provide the proof for completeness. 

Proposition A. 4. Let the flux f be uniformly convex, i.e. f" > c > 0, and for u G 

L°°(R) let u G L°°([0, +oo[xR) he the solution provided by Theorem \AM Then u G 

L-(]0,+oo[;SMoc(K)). 

Proof. We decompose the measure d^uit, ■) as 

dxu{t, ■) = [dxu{t, ■)Y - [dxu{t, ■)]" . 
Then, for every i? > we have 

\dxu{t,-)\{[-R,R]) = [dxu{t,-)Y{[-R,R])+[dxu{t,-)y{[-R,R]) 
= 2[dxu{t,-)Y{[-R,R])-[dxu{t,-)\{[-R,R]) 

< ^ + \u{R+)-ui-R-)\ < l^ + 2||n|U. 
ct ct 

□ 

The following results regarding the strong continuity of the entropy admissible solution 
has been proven by Chen and Rascle [T6] . 

Theorem A. 5. Assume that f"{u) > c > for any m G M and let u G L°°(]R). Let 
u G L°°([0, +oo[x]R) be a solution in the sense of distributions of flA.6p with n(0, ■) = u, and 
assume that u is entropy admissible. Then, u G C° ([0, +oo[; -^[^^(IR) — s), and in particular 
u is the unique solution provided by Theorem \A.2[ 

Remark A. 6. Proposition IA.4I and Theorem IA.5I can be easily extended to the case of a 
uniformly concave flux, that is /" < — c < 0. This can be seen by setting v{t, x) = u{t, —x) 
and observing that v is an entropy admissible solution of the equation dtv + dx \{,—f){v)\ = 0, 
which has a uniformly convex flux. 

Remark A.7. If /' is globally bounded, then Theorems IA.2I and IA.3I and Proposition IA.4I can 
be extended to the case of a locally bounded initial datum u G Lj^^(]R), the only difference in 
the results obtained being that now the solution u belongs to L'^^(^[0, +oo[xM). Indeed the 
constant L in (lA.SP and (]A.9P actually depends on 1 1 w 1 1 and llt'llioo only through ||/'( 
and (see (1A.7P ). Being these quantities finite, this ensures a global bound on the 
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speed of propagation. In addition, observe that the notions of solution in the sense of 
distributions and of entropy admissible solution are both local. Hence, by considering a 
suitable truncation of the initial datum and exploiting the finite propagation speed, we can 
easily complete the argument. 
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